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by 

J.  Dunham 


SUMMARY 

A  historical  review  is  presented  of  the  influence  of  computational  fluid 
dynamics  on  the  developsient  of  UR  compressors  and  turbines.  The  ability  to 
predict  pressure  distributions  has  led  to  increasingly  successful  attempts  to 
tailor  aerofoil  shapes  in  such  a  way  as  to  optimise  performance.  Once  proven, 
new  computer  programs  have  rapidly  been  put  to  use  by  Rolls-Royce  and  by  same  non- 
aeronauticsl  firms.  Examples  are  given  of  improvements  achieved  hy  their 
application. 


The  present  state  of  the  art  is  assessed,  and  the  prospects  for  future  com¬ 
putational  developments  are  discussed.  The  need  for  detailed  experimental  test 
cases  is  emphasized.  The  economical  representation  of  viscous  effects  remains  a 


key  difficulty,  especially  when  heat  transfer  predictions  are  needed.  (  ^ 
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On*  of  tlM  aayaadlag  imi  af  scientific 
research  In  th*  logo's  haa  boon  soayntatloaal 
fluid  dynaalcs.  CTO  ptaMili  a  aaal  stlMlatlag 
Intallaataal  ehallaaga  to  tha  Mthaaotlcally- 
glftad,  aad  affara  a  any  far  tha  fluid  dyaaalelst 
ta  daaalaa  Ui  art  daaylt*  tha  lMreaalag  coat  aad 
tftatala  af  high  syaod  nndynaalc  axy*rVa**ti ■ 
■at  daaa  C71  nelly  half  tha  daolgaar  pro dean 
hattar  saglaaarlag  ftahita,  or  radaat  davalaaaaat 
ttaa  aad  caatf 

Tha  fttfiii  af  thla  pafar  la  ta  gin  a  claar 
afrtraatlaa  rafty  ta  that  ftaatlaa,  ha aad  *aaa 
raotawtag  tha  Mat  ary  af  aaraaa(laa  tarts** shlasry 
raaaareh  la  tha  R  aaar  tha  paat  19  yuan.  It  la 
traa  that  other  flalda  aaa  ha  aaaa  la  ahlch  CTD 
haa  aat  yet  daaaaatrahly  halyad.  A  atady  af  tha 
yartlaalar  alraaaataaaaa  af  thla  hlataty  laada  ta 
aaaa  ■■aural  eaatlaalaaa  ahaat  taadltlaaa  aaadad 
far  tha  aaaaaaafal  daaalayaaat  aad  cnylaltatU*  af 
ere. 


Tha  yayar  la  written  fro*  tha  yolat  af  elan 
af  a  “aaar"  af  tarhaaaehlaary  CTO  rathar  than  a 
dawalayar  af  ltt  aa  attaayt  ta  aada  to  raalaw  tha 
analytical  aad  eaayatatloaal  aahaaaa  no ad,  or  ta 
auk*  Jadgsaaato  ahaat  thalr  talatlaa  aclaatlflc 
■art  to. 

It  la  loan rt aat  ta  note  that  la  cltla«  thla 
BK  caaa  history,  It  la  aat  auggaacad  that  tha  OK 
haa  aalaaaly  achlaaad  CTD  aaccaaa.  ilallar  report 
al*ht  wall  ha  wrlttaa  ahaat  aerroayaadlag  axyarl- 
aaea  la  ether  ceaatrtee.  Crahaa  ddaaaayh,  aad 
■ohllk'1'  her*  already  dona  aa  ror  the  DU. 

The  yayar  haylaa  by  outlining  tha  broad 
yattore  af  euelwtloa  of  aaraaaglaa  tarhaaaehlaary. 
A  historical  rorlau  of  tha  daaalayaaat  of  turbo- 
aachlaary  ere  aat  hade  fellow*.  Th*  re  la  thaaa 
aathada  haa*  ylayad  la  tayroalag  aaraaaglaai  la 
thaa  daacrlhad,  with  aa  My  laa;  further  aaaaylaa 
are  taka*  flea  tarhaaaehlaary  la  ayylleatlaaa 
ether  thaa  aara  aaglaaa .  After  eaaaaata  aa  thla 
hlatary,  tha  fatar*  aaada  aad  yaaalhllltlM  at 
tarhaaaehlaary  CTD  are  coaaldarad. 


It  ta  flrat  aacaaaary  ta  ahaaraa  tha  aatard.  ( 
af  aara  aaglaa  daaalayaaat,  Mich  la  slgnlfloaltly  t 
dlffaraat  frua  that  af  aircraft  or  of  yaar/ 
gaaaratlaa  ylaat.  Aa  aaglaa  la  flrat  doaTgaCd 
ualag  tha  heat  tachaalagy  aoatlahla  at  tha  tlaa, 
aad  aooa  yratocyyoa  are  aada  and  taatad.  Deeytte 
tha  haat  efforts  of  tha  daeigaare,  thaaa  ymew- 
tyyaa  exhibit  daf acta  la  th*  fan  af  yraaataru 
fallaraa  or  laadagaata  yarfensaae*.  Far  savaral 
year*,  tataaaloa  daaalayaaat  thaa  follow*, 
eerrwctlog  daf  acta  aatll  a  faraal  tyya  ayyrooal 
tat  ta  yaaaad  aad  aaglaa*  aru  put  late  aarwle*. 
Th*  daaalayaaat  yroaaaa  daaa  aat  and  than, 
frograealoa  layrwanaaata  ta  hath  durability  aad 
yarfsraaars  caatlao*  to  ha  triad  aat  on  deaaloy- 
ant  aaglaea,  aad.  If  ayyrayrlat*.  laearyoratad 
tat*  yraduetloa.  Oaually,  th*  Might  af  aa  air¬ 
craft  las  re  eon  la  auccaaatr*  oarslana  a*  sirs 
yaaaaagars  an  prwoldad  far  or  nan  uaayoas 
carried;  thla  ragaira*  aara  passer  f raw  th* 
aaglaaa.  Th*  icaenlcs  af  hath  elatl  aad  allltary 
alnnft  ayantloM  jMttfy  far  aan  rhaagas  la 
aaglaa  hladlag  thaa  aaald  soar  ha  Ratified  ta, 
•ay,  tha  daalga  af  ••  alnnft  wing.  la  th*  yarn 
g aaar ac laa  industry,  tha  gaaatlclaa  af  ataua 
turhtaaa  ar  yaa  turhlaaa  aald  an  aat  (officiant 
t*  Ratify  a  danlayaoM  yngnaas  af  the  aoro- 
•aglM  tyya.  It  la  claar,  thanfan,  chat  aarw 
aaglMt  haa*  a*ah  aan  fragaaat  syysrtualtlaa  to 
latradaaa  aa*  taahaalagy  thaa  aaay  achat  fluid 
dyaaada  dawlaaa. 

Leaking  back  at  tha  hlatary  of  oanaaglM 
tarhaaaahlaary  ta  tha  OK,  it  la  yaMlhla  to 
Idaatlfy  a  asll  auabar  of  ajar  ados  or  is  in 
aaayraaaar  taabaalogy  and  tarhln  tachaalagy 
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which,  In  cooblnatien  with  mv  iccnlaU  end 
ntnufteturing  procwcwwa ,  hen  cublad  n|lu 
aiaadarda  to  wove  thud.  In  hitmn  thaaa  we  apt, 
the  prog rvttlv*  refloenent  of  individual  deolgn# 
haa  tap  raved  chan  to  the  Unit  o<  tha  currant 
technology.  Developaant  la  tha  procaaa  of  clia Ir¬ 
ina  up  tha  aapaptocaa;  raaaarch  la  tha  procaaa  of 
aovlof  tha  ttpwptottc  thanaalvaa* 


The  broad  pattam  of  avolutloe  of  axial 
cowpreeear  technology  in  tha  OX  la  illuacratad  in 
Fig  1.  hath  tha  raaaarch  dtraaa  (on  tha  laft)  and 
tha  davalopaaat  otraaa  (aa  tlUL  right)  ocarted  t ran 
tha  aarlv  work  of  Griffith'2',  Conataot'J',  and 
huwull'21  in  tha  ltlO’a  and  IMO'a.  Tha  aajor 
achlaaaaaat  of  tha  ItSO’a  aaa  tha  Up  roved  uader- 
ataadlng  of  ataga  Batching,  a tall,  rotating  a tall, 
aad  blada  vibration  aa  which  awlti-epool  aaglnaa 
warn  baoad.  Tha  naxt  aajor  achiavaaaat  aaa  tha 
avolutlaa  la  tha  lata  IMO'a  at  tha  atraaallaa 
cwrvntwra  not  hod  far  aaalyalag  anrldlooel  flaw  (of 
which  worn  latar)  ahich  allawad  daaigna  to  ha 
praparad  along  a t ratal  1  aaa,  aad  aa  aaablad  tran- 
aoalc  faaa  to  ha  daatgaad  aad  davwlopad  to  far 
higher  ataadarda  of  efficiency  than  pravtoaaly 
achlavad,  ho  aajar  aaa  aaradywaalc  advaacaa  caaa 
aftar  than  aatll  tha  IMO'a  whan  tha  idaaa  of  "aad 
haada"  (local  profile  aodlflcatlaaa  aaar  aad 
walla)  aad  aoparcrlclcal  profliaa  wara  flrat 
appliad  aad  whan  it  baraaa  peaalbla  to  lndwda 
vlacoua  affacta  in  blada-to-blada  calculat loaa . 

The  origin  of  OK  centrifugal  coapraaaora 
(Fig  2)  Una  In  tha  awparchargara  at  aircraft 
platan  auglneo  batata  a  tha  wara.  Thaaa  wara 
davalopad  bp  Whittla'a  taaa  to  a  atata  of  relle- 
bllltp  aad,  for  tha  tine,  high  efficiency'1'.  A 
aajar  atop  forward  waa  aada  poealbla  la  tha  IMO'a 
bp  tha  avolutlaa  of  nee  bod  a  far  tha  aaalpala  at 
lnvlecid.  aerldteaal  flan  aad  thaaa  famed  tha 
baale  of  coaaidarabla  advaacaa,  eepaclallp  la  the 
OS*  and  Canada.  inapt-back  bladaa  ware  Intro¬ 
duced. 

Turblnaa  far  aarlp  Jat  aaglnaa  (Fig  3)  wara 
dealgaad  waing  a  taaa  turbine  practlaa,  but  adapted 
far  vertex  flaw  bp  Whit  tin  aad  ethara.  Tha  flrat 
aajar  raaaarch  advance  amaa  free  *lalap  aad 
Mathleeoa'a  part oruaaca  pradlctiau  Bached<*5, 
bacauaa  it  aaablad  aptlawta  walacltp  trlaaglaa  to 
be  talactad  far  tha  tpaclf lad  daep.  Tha  naxt 
elgatficaac  advaaaa  earn  uhaa  aathoda  at  calculat¬ 
ing  tha  velocity  dlatrtbatlaa  a  round  the  blada 
aurfacat  were  auployod.  Thao  latar  cam  tha 
■  t  meal  In#  curvature  through- flaw  uathed,  particu¬ 
larly  helpful  for  turblnaa  with  highly  flared 
anawlue  liana,  haw  advaacaa  am  nan  aaargiag  la 
tha  IMO'a,  In  tha  fora  af  aad  wall  chap  tag  ar 
blada  etaaklag  rhaagaa  with  tha  abject  af  reducing 
aacaadary  laaaaa. 

Looking  back,  neap  at  tha  aajar  advaacaa  la 
t urbane ahlaarp  technology  aver  tha  loot  thirty 
paara  have  occurred  whan  daalgaero  wara  able  to 
aatlaata  praaawra  dial rlbwtleaa  an  wetted 
aurfacat,  aad  aa  awpreoa  gwaatltatlvely  tha 
Intuitive  aadantaadlag  at  turbaaachlaary  flam 
which  rvaaarth  wavfcara  davalopad  flrat.  The 
aathoda  far  delag  aa  will  now  ha  r art awed. 


FXGUKS  1 .  Axial  Conpraaaor  Hiatory 


ntniaacM  DcvnusatagaT 


riGUag  2.  Centrifugal  Caagiraaaor  hiatory 


riCOSS  3.  Axial  Turbina  History 


Tarhsaachiaery  haa  always  beea  gaatgaal  aa  a 
~quaal-thraa-4laaaaloaat  naiiy  (laa*  heels 
(Fig  4).  IW  geelga  process  stares  by  abaaa lag  a 
~thrsagb-(lew*  patters  aa  a  aarMiaaal  glass  -  a 
crooe-asctlaa  of  tha  aaglaa  laclaglag  tba  aala  ml 
ratatlaa.  tbs  tbrsagb  flaa  la  aaaealaeg  aa  a 
clrcaalaraaclal  asarags  ml  tba  raal  (laa,  at  aa  a 
typical  a  erase  sarfaca  pass  lag  bstassa  tba  blaSaa. 
Sacaagly,  tba  Saalga  al  tba  blaSa  profiles  (araaa- 
aactlaas)  la  aaaalSatag  la  a  bl  sSa-ts-blaSa- 
atraaa  aarfaaa,  if  lag  tba  (laa  ta  ba  ataagy  la 
tlaa.  tbta  lag  1  laa  that  tba  saataaMaaaa  ml  tba 
(law  least  ag  prsstasa  raws  la  aa44aaly  t  las- 
aearagag.  Tbta  gaaaral  ayyrsacb  aaa  Itrat 
(oraally  jsstiflag  by  He<*>  la  IWI. 

Tba  thrssgb  (laa  calc  slat  laa  drat  gaaarataa 
tba  aalayaaatrle  atraaa  sar(aaa  far  tba  blsgr-ta- 
blags  yragraa.  That  la  tars  gaaarataa  tba  tetal 
yrmssra  laaa  aag  tlaa  gtraatlaa  ragalrag  ta r  tba 
tbraagb-daw  yragraa.  rtaally,  attar  Itaratlaa  It 
aaagag,  tba  tbreagt-llaw  yragraa  salaalataa  tba 
esatall  abaraatarlatlca  at  tba  tarbatablan  (law, 
y  rasas rs  ratls,  aag  a((lclaaay.  Tba  rale  at  tba 
blaga-ta-blaga  yragraa  la  ta  aayyly  leaal  (las 
•agio  aag  yraaaara  lass,  aag  ta  yraglat  blags 
aer(eaa  yraaaara  giatrlbatlaai  aaaa  yragraas  a Isa 
yraglat  baaatary  layara.  Tba  rale  at  tba  tbraagb- 
daw  yragraa  Is  ta  yraglat  raglal  (law  glatrlbw- 
tlaa  aag  aaarall  yarfsraaasa. 

Vatll  ralatlsaly  rsaastly,  tsrbaaaablaary 
analysts  taalg  aaly  ba  gaaa  aa  tbla  gaeel-tbree- 


neuu  4.  Tha  Qoasl-Threa-OUaanslonal 
Approach 


glaaastaaal  staagy  (law  baaia,  bat  mw  asny  tally 
tbtaa-glaaaalaaal  aag  aaaa  aastaagy  (law  caaywta- 
tlaaa  are  yasatbla.  It  la  therefore  caasaalaat  to 
review  tba  history  sf  tarbaaachlaery  CFO  eager  tba 
eab-beattags  ml  bl«ga-tw-blage  aacbaga;  tbraagb- 
daw  watbaga;  gsaet-M  aathoga;  aag  felly- ID 
aat hags.  For  tba  yaryaas  o(  tbla  yayar,  a  CFS 
aatbag  la  ga(l sag  aa  a  aalatlaa  ml  tba  lsslaclg  or 
slaaaaa  agaatlsaa  sf  (led  aotlsa  la  two’  ar  throe 
glaaasloaa  (as  glstlaat  (raw  oaa  glaaaaloa). 
latagral  baaagary  layar  aarbaga  sag  bast  traao(ar 
Bethels  are  easlagag  ta  baay  tba  laagtb  ml  tba 
g  1  saws s laa  wltbla  baaaga. 


Tba  yraaaaa  aa  (ar  gascrlbag  la  tba  calcsla- 
tlsa  at  tbs  (law-dalg  wltbla  a  tsrbsas china  ml 
aysalflsg  ysaastry.  It  ta  alas  gaatrsblo  ta  bare 
CFD  gaalga  yragraas,  ablab  will  caayate  tba  shays 
ml  blaglag  ragalrag  to  aehlaoa  a  syaatdag 
yartaraaaaa  ta  aa  oytiaaa  way.  Design  yragraa 
are  aaaalgarag  eager  a  dasl  ssb-baaglag. 


Thors  tea  baft  a  acaagy  awaaaailn  a t  haabe(*~*>, 
yayara'1®"**' ,  aag  Car?araaaaa  aag  tartars 
»arlas<,0_“>  gasatag  aatlraly  ta  tsrbsas shlaary 
Cft.  Tba  rasisr  la  rsferreg  to  tbaa  (ar  gatallag 


aayaaltlaaa  at  tba  agaatlsaa  aaag,  tbs  aalatlaa 
aeahagt  sgaycag,  sag  tbs  jasctdaaclsa  at  tba 
rasalta  by  saayarlssa  slth  aysaldc  aayargaaats. 
Tba  yaryaas  al  tba  yraaaat  yayar  ta  ta  steal  hash 


free  tba  Sat  cl's  aag  rastaa  tba  historical  aaa  aag 
sllaatlaaaaaa  at  tba  watbaga. 
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ilade-to-blsde  — thods 

TEa  urunt  calculations  of  flow  through  a 
cascade  of  cambered  blades  of  flalea  thickness 
wars  by  Bow«ll<23),  and  Ha r chant  aad  Collar^2*), 
both  la  1941.  They  sol  rod  eh«  lavlscld  ineow- 
prssalbls  steady  plane  flow  equations  (which 
reduce  to  Laplace's  equation)  by  eonforaal  trans¬ 
formation.  Many  others  followed  (see  p268  of 
ref  7)  Including  the  extension  to  linearised 
coupreoelble  flow.  The  alternative  c leas leal 
■at  hod  of  solving  Laplace's  squat  loo  is  the 
singularity  nathod,  aad  assy  ouch  aathods  warn 
evolved  because  they  eould  analyse  arbitrary 
profiles.  The  first  method  adopted  In  the  WL  for 
practical  uoe  was  that  of  Nnrtensen'23',  adapted 
for  eeapresslble  flow  by  Price'**'. 

In  the  1940's,  wore  versatile  ways  of  solving 
compressible  flow  equations  were  developed,  as  the 
growing  power  of  computers  aade  such  aathods 
practicable.  They  were  generally  elaborated  to 
solve  the  sebeoulc  eeapresslble  flow  around 
cascade  blade  profiles  oa  a  streea  surface  of 
arbitrarily-varying  thickness  and  radius,  as 
required  for  the  quasi- 30  approach.  Four  types  of 
method  evolved!  finite  difference,  finite  differ¬ 
ence  using  matrix  Inversion,  streamline  curvature, 
end  finite  elaaeat. 

ly  subsequently  employing  say  convenient 
boundary  layer  as thud,  predict leas  of  total 
pressure  loss  eould  be  attempted.  The  varying 
streantube  thickness  aad  radius  eould  easily  be 
catered  for  in  ea  integral  boundary  layer 
nethed'27'.  Outlet  angle  prediction  proved 
troublesome,  because  the  trailing  edge  of  a 
practical  terbemechlae  blade  is  not  sharp  eoo^h 
to  apply  the  Kutta  condition  simply.  Instead,  a 
convergence  ef  surface  pressures  towards  a  cannon 
trailing  edge  presence  mss  adopted. 

There  use  mush  discussion  at  tbs  time  ef  the 
relative  ce^utetiemsl  aad  practical  merits  of  the 
four  types  of  method  listed.  The  amount  of  detail 
needed  at  the  leading  and  trailing  edges,  especi¬ 
ally  at  off-design  Incidence,  was  highly  relevant. 
In  the  UK,  the  streamline  curvature  method  use  the 
only  sms  to  emerge  at  this  time  as  a  regal  sr 
design  tool.  It  is  probably  true  that  this  urns 
doe  not  se  much  to  the  relative  scientific  merits 
of  the  fear  methods  as  to  the  fact  that  tolls- 
Reyce  had  invested  a  large  amount  ef  effort  within 
the  Company  into  developing  st  re  online  curvature 
programs  into  practice!  tools  capable  of  being 
applied  by  specialists  other  than  the  program 
originator. 

The  singularity  and  matrix  blade- to-b  lade 
methods  have  also  bepn  applied  to  radial  turbo- 
machines,  by  lailly'**'  and  ©onion'**'  respect¬ 
ively. 

Although  the  methods  se  far  described  sen 
generate  solutions  with  supersonic  patches,  modern 
terboanchinoo  centals  check  moves  and  supersonic 
regions.  It  nos  therefore  a  grant  step  forward 
when  goatee  developed  practical  ^  time-nnreblag 
methods  for  turbenaehlnos'30' ,  following 
NecCerneck's  approach,  which  ware  cepeble  ef 
computing  flew  fields  with  embedded  shock  waves, 
he  used  a  finite  volume  formula  time,  which  he  has 
since  substantially  improved'31'.  Ms  aathods  are 
new  used  by  many  firms  throughout  the  world. 


Still,  the  flew  through  a  compressor  cascade 
with  supersonic  inflow  could  not  be  realistically 
noddled.  The  leading  edge  shock  impinges  on  the 
suction  surface  end  either  separates  the  boundary 
layer  or  brings  it  elose  to  separation.  The 
natural  way  of  modelling  the  boundary  layer  is  to 
compute  the  boundary  layer  displacement  thickness, 
add  it  to  the  blade  natal  thickness,  recalculate 
the  velocity  distribution,  and  Iterate  to  convert 
once.  Unfortunately  any  attempt  to  do  so  is  found 
to  be  unstable  with  a  separating  or  nearly  separa¬ 
ting  boundary  layer.  This  difficulty  woo  at  last 
overcome  by  Calvert'3*'  la  I9t2,  by  inverting  the 
order  ef  iteration  along  the  suction  surface  after 
the  shock  wave.  Is  obtained  converged  solutions 
tfi  cases  with  separated  or  attached  boundary 
layers  which  agreed  well  with  experiments'33'. 

None  of  the  nethode  so  far  described  ere 
mined  et  computing  unsteady  flow.  VhlteheadC3*) 
bee  developed  e  finite  element  method  for  this 
purpose,  which  Has  been  used  ee  e  design  'end 
analysis  method  for  etoedy  flow,  ee  well  ee  form¬ 
ing  the  heels  for  e  flutter  predlctlom  system. 

These  methods  of  Don tom,  Calvert,  end  White¬ 
head  have  been  adopted  by  Aolle-toyce,  extended 
where  aeceseary,  end  applied  regularly  to  analyaa 
and  deelgn  compressors. 

Attentloe  bse  slso  turned  in  the  19§0*e  to 
solving  the  viscous  squat ions ,  now  thst  super- 
conputtrs  are  norm  widely  available.  The  problem 
le  how  to  model  the  turbulence  economically  and 
yet  sufficiently  realistically,  for  the  turbo- 
nechlnery  application.  Dawns'391  adopted  a  simp Is 
nixing  lsngth  model  with  some  success. 
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nearly  slaultsaeously  by  nod 
Novak'403  la  the  USA  and  by  Hetherlngtoa'413  and 
ftingroee  in  the  UK.  Ia  lea  original  fora,  calcu¬ 
lation  planes  wara  located  at  every  trailing  adga 
plana,  w  that  local  outlet  angles  aad  total 
preaeura  loaaaa  could  be  lnaartad  froa  caacada 
corralatioaa.  The  Hetherlngtoa  aathod  was  adopted 
by  tolla-toyce  and  haa  bean  widely  uaad  aver  a Inca 
for  the  design  and  analytic  of  both  compressors 
aad  turbines. 

The  addition  of  calculating  planes  within 
roue  caused  ouch  discussion  In  the  early  1970* a, 
when  it  harsmo  dear  that  the  radial  coopoaents  of 
blade  force  (not  accounted  for  when  only  trailing 
adga  plaaaa  are  used)  could  lafluaaca  the  flow 
considerably.  Smith'4*3,  using  a  aatrix  as t hod 
with  a  grid  within  the  blade  rows,  demonstrated  a 
auch  batter  prediction  of  the  radial  variation  of 
•tatle  pressure  in  a  turbine  (which  would  affect 
coollna  airflow  predictions).  Nora  racantly. 
Cinder!433  haa  shown  the  laportanca  of  tho  antra 
plaaaa  la  a  streamline  curvature  prograo  applied 
to  a  transonic  fan. 

Through-flow  methods  have  also  beta  developed 
for  radial  turbooachlnee.  la  1967,  Hood  aad 
Marlow!443  applied  the  atraaallae  curvature  aathod 
to  a  pump  legalist.  At  the  National  Use  Turbine 
Establishment,  en  'lonelier  Cooputer  Design 
Package**  wee  evolved'433  which  combined  a  matrix 
through-flaw  calculation  (with  an  approxlaate 
blade- to- blade  assumption)  with  stressing  and 
numerically-controlling  machining  elements.  At 
the  Nstlooal  engineering  Laboratory  sn  Incompress¬ 
ible  flow  pump  duelgn  package  was  produced'443, 
alao  os  lag  a  matrix  method,  aad  asking  use  of 
lallly'a  singularity  aathod  for  blada  section 
design.  Coulee'47'  haa  also  developed  the  aatrix 
method  for  centrifugal  coapraaaors. 

One  of  the  features  of  multistage  coopresaor 
porforaaoee  which  levlacld  through  flow  methods 
cannot  predict,  and  which  through! low  aethoda  with 
annulus  wall  boundary  layer  allowances  also  fall 
to  predict,  la  the  “repenting  stage"  phenomenon 
identified  by  L  ■  Smith'4®3.  The  redial  distribu¬ 
tion  of  axial  velocity  settles  down  after  a  few 
stages  to  a  fairly  constant  pattern.  Ia  1981 
Adkins  aad  Smith'43 J  proposed  an  explanation  knead 
essentially  en  eacoodary  flow  phenomena,  end 
introduced  semi-empirical  teems  late  their  at  rear 
lime  curvature  program  to  calculate  the  secondary 
flow  affects.  Their  results  were  consistent  with 
experimental  observations.  Nora  recently, 
Csllieors'^J  boa  proposed  aa  alternative  explana¬ 
tion  booed  an  nixing  theory.  This  provides  aa 
avaa  aore  convincing  explanation  In  the  particular 
caaaa  tasted,  though  aa  la  all  alxlag  theories  tbs 
prediction  of  affective  turbulent  eddy  viscosity 
la  problematical.  The  full  lmpllcatloas  of  mixing 
on  turbemechlaery  dealga  have  yet  to  be cams  clour. 

Until  larger  aad  faster  computers  were  devel¬ 
oped,  it  wee  essential  to  avoid  having  to  solve 
fully  three-dimensional  agnations  with  the  proper 
boundary  conditions.  One  approach  wee  fundament¬ 
ally  aa  actuator  disc  approach,  la  that  fully 
three-dimensional  solutions  were  obtained  on Jv  in 
unbleded  ducts  (ftawtherae' 51  3 ,  Dunham'5*3, 
Lewi  s'53').  This  approach  has  net  lad  to  a 
practical  werklag  method. 


The  universally  adopted  method  was  first 
suggested  in  1951  by  who  derived  the  basic 

equations.  It  involves  Iterative  solution  of 
blede-to-bladu  flow  end  through-flow,  as  described 
earlier.  This  process  converges  in  only  a  few 
1  care  dons,  to  principle,  it  does  net  allow  for 
the  effects  of  streeawlee  vortlcity,  which  die- 
torts  the  assumed  blade- to- blade  stress  surf sees. 
Coalas' 34 '  baa  proposed  extra  tarns  for  atream- 
wlse  vortlelty  of facta  In  centrifugal  compressors. 
However,  the  basic  method  le  used  for  design 
purposes,  aad  haa  bean  limited  primarily  by  the 
quality  of  the  blede-to-blade  end  through-flow 
methods  incorporated  into  it. 

fully  three-dimensional  aethoda 

sines  yJUft  and  Httharlngton'333  end  Oliver 
and  Sparis*543  first  proposed  fully  three- 
disens looal  aethoda  lo  1970/1,  at  least  sixty 
papers  have  bean  published,  moat  of  then  based  on 
tho  KacCormack  tine-marching  approach.  The  main 
contributors  of  lavlacld  echenea  have  bean 
Denton'573,  Thompkina'5*3,  »©eman'533  (who  all 
first  published  in  1976),  aad  Wreck'403  (1980). 
The  faatnra  of  all  these  aethoda  la  that  the 
convergence  in  tins  la  alow  -  needing  of  the  order 
of  1000  etepo  -  and  this  haa  to  bo  controlled  by 
select  log  the  rljht  grid  ahape,  end  the  beet  tlae 
step,  ead  the  beet  relaxation  factors.  The  Denton 
finite  volume  fomlatlon  haa  bean  laproved  and 
extended!**3  over  the  years  since  1976.  This 
numerical  schema  conserves  mass  flow  exactly  but 
does  not  automatically  conserve  etagnetlon 
pressure.  Thompklns  uses  e  finite  difference 
scheme  and  Hlrech  a  finite  element  schema. 

One  of  the  first  attempts  at  a  fully  three- 
dimensional  viscous  calculation  appaere  to  be  that 
of  Carrick  (1975)'423  who  added  simple  viscous 
tsras  to  Stuart  and  Hatherlngton" s  aathod'333. 
This  wus  not  pursued,  however,  and  tha  major 
contributions  hava  been  made  by  tfalltt  (froa 
1975)'633,  Dodga  (froa  1977) C443  and  tha  Moores 
(froa  1979)<45>. 


The  tfalltt  and  Dodge  methods  aro  essentially 
successive  approximations.  Tha  equations  srs 
written  with  lavlsdd  tsras  on  tha  laft  hand  aid# 
and  viscous  tana  on  the  right  head  aide.  Start¬ 
ing  froa  an  lnvleeld  solution,  tho  viscous  tsras 
can  then  be  calculated,  end  further  solutions 
obtained  successively.  Beth  axial  aad  centrifugal 
flow  fields  have  been  calculated  by  theea  aethoda. 
The  Neoree  use  a  finite  difference  scheme,  vhleh 
le  coaflaed  to  eubeonle  flow  la  principle.  An 
laprovad  version  has  recently  bona  devised'4®3. 

Keceatly,  Denton^413  aad  Dewee^333  havn 
proposed  separate  schemes  for  extending  the  tiae- 
aarehlng  approach  to  viscous  floae.  Dewae  nose  a 
alxlag  length  nodal  ae  ia  hie  twa-dlmeaaloaal 
cede,  aad  Deaton  haa  tried  avoiding  modelling  mddy 
viscosity  by  using  aapirical  blade  force  teems. 
The  eeeeacUl  difficulty  with  ell  the  viscous 
aethoda  le,  of  course,  accommodating  a  suffici¬ 
ently  fine  grid.  To  resolve  the  flea  sear  wells 
it  is  essential  to  have  sufficient  grid  Hues 
actually  within  the  boundary  la yurt  this  enquires 
•  huge  computer  store. 


Doslao  as  s  dlrect  PfPCSM 
No  formal  solution  has 
design  problem  das  it  Is 
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general  tenmt  conpute  Che  optlnun  turbouechine 
for  a  specified  duty.  The  number  of  Independent 
variables  le  effectively  Infinite,  and  the  maber 
of  design  object ivee  le  conelderable  and  their 
relative  importance  unquantlfled  (In  general)* 
These  design  objectives  can  be  stated  as  follows* 

1.  smchanleal  Integrity 

This  sets  Units  to  leading  and  trailing 
edge  thickness,  and  thickness  In  general 
(especially  for  cooled  blading).  Xt  also 
requires  freed on  froe  significant  vibra¬ 
tion  of  any  kind  over  the  running  range. 

2.  efficiency  over  the  running  range. 

3.  surge  margin  (for  a  compressor)  and  low 
susceptibility  to  inlet  flow  distortion, 
over  the  running  range. 

A,  low  else,  weight  and  cost. 

5,  minima  design  and  development  cost, 

When  undertaking  a  design,  proposals  must  be 
assessed  against  aach  of  these  objectives  to  as  to 

roach  a  final  compromise;  this  Is  still  e  human 
Judgement. 

tt  le  necessary  to  start  by  making  an  arbi¬ 
trary  trial  choice  of  most  of  tbs  Independent 
variable*.  In  practice  this  la  dons  by  selecting 
s  through-flow  pattern  conforming  to  optimum 
volume  of  etege  loading  parameter*  known  from 
previous  experience;  that  le,  the  "optimum*  veloc¬ 
ity  triangles  are  selected.  Any  eoaputar-besed 
optimisation  Is  undertaken  In  the  blede-to-blede 
calculation ,  where  the  designer  ac tempts  to  find 
the  optimum  blade  profile  for  the  required  inlet 
and  outlet  flow  angles  sod  velocities.  Sven  this 
is  not  easy. 

Host  approaches  Involve  prescribing  e  "good* 
velocity  distribution  (?VD)  end  computing  the 
blade  profile  required  to  generate  it. 
Llghthllil*7'  proposed  e  FVD  method  in  1949,  a* in* 
e  conformal  transformation.  In  1952,  Stanlts'*®' 
Introduced  e  linearised  compressible  FVD  method, 
which  wee  developed  by  Faynet**)  sod  applied 
within  Rolle-toyce  to  design  turbine  blade  shapes. 
Hurugesen  end  Ralllyj'j)  wrote  e  design  version  of 
Her  tensin' a  methods' end  lewis*71'  also  wrote 
s  design  method  using  distributed  singularities. 

The  fundament  el  difficulty  in  using  the  FVD 
approach  is  chat  the  selection  of  the  pressure 
distribution  requires  experience  (end  trlal-aad- 
error)  If  an  unacceptably  thin  aerofoil  (or  evou 
one  of  "negative  thickness",  since  the  Stanlts 
method  actually  designs  s  peasage)  la  to  bo 
avoided.  Wilkinson'72'  adopted  an  Interesting 
schema  to  avoid  this  difficulty;  hie  method 
designed  the  section  surface  (with  the  more 
critical  velocity  distribution)  to  satisfy  e 
prescribed  velocity  distribution  hue  then  pre¬ 
set  l bed  the  aerofoil  thickness  end  calculated  the 
result lag  pressure  surface  velocity  distribution. 
A  streamline  curvature  method  on  those  lines  has 
been  developed  within  Rolls-Royce,  end  used  to 
deolgu  turbines. 

The  levers  ton  of  a  singularity  method  into 
FVD  form  involves  solving  the  seme  equations  with 


different  Variables  unknown.  The  inversion  of 
more  modem  analysis  methods  (using  grids)  is  mors 
complicated  and  involves  iteration  of  the  analysis 
program.  Paige' 7  27  has  written  e  FVD  finite 
volume  method  using  a  hill-climbing  scheme.  Cedar 
end  Stow'7*'  hevo  written  e  FVD  finite  element 
method  based  on  Ref  (34),  using  e  local  transpira¬ 
tion  nodal  to  avoid  changing  the  grid. 

Recognising  that  trial-and-error  Is  necessary 
In  applying  even  a  FVD  schema,  a  workable  alterna¬ 
tive  la  to  provide  a  flexible  shape  description  In 
the  form  of  one  or  moro  algebraic  express  lone 
Involving  arbitrary  parameters.  Than  the  choice 
of  shape  is  controlled  by  the  choice  of  a  fixed 
number  of  paramotors,  say  9.  Starting  from  guide- 
linos  established  by  experience,  the  parameters 
arc  varied  until  a  satisfactory  pressure  distribu¬ 
tion  and  machanlcal  shape  are  obtained.  A  method 
of  this  evpe  for  turbine  design  wss  proposed  by 
Dunham*75'  end  this  approach  Is  also  used  by  RAE 
for  transonic  fans. 

A  more  radical  approach  was  proposed  by  la 
Foil'7*'  who  took  the  further  step  of  prescribing 
the  desired  boundary  layer  development  end  hence 
working  beck  to  the  profile  shape.  This  net hod 
would  presumably  encounter  the  seme  difficulties 
as  FVD  lm  arriving  at  a  machanlcal  ly  accept  able 
•haps;  it  has  apparently  not  been  adopted  by  any 
manufacturer. 

4  Apeli cat loss  o  f  CFP  to  Improvlna 
Turbmochinerv 

lm  this  sect  lorn,  the  historical  role  of  CFD 
In  improving  turbomachinery  la  examined,  end 
Illustrated  by  some  specific  examples  taken  from 
aeronautical  and  Industrial  applications. 

Axial  Compressors 

Early  axial  compressors  wars  designed  using 
*atandard"  profile  shapes  evolved  from  systematic 
cascade  testing,  end  assuming  simple  radial  equil¬ 
ibrium  conditions.  Although  Ineonprueotblo  blade- 
to-blade  design  end  analysis  methods  were  avail¬ 
able  lm  the  1940's,  and  later  compressible  ones, 
they  were  never  used  to  design  engine  blades, 
because  of  their  inability  to  predict  the  key 
features  of  outlot  flow  angle  end  total  pressure 
loss.  The  compressors  used  In  engines  designed  In 
the  1940's,  1950's,  sod  early  1990's  were  devel¬ 
oped  to  acceptable  performance  end  reliability  by 
naans  of  long  expensive  test  programmes.  Many 
cascade  tests  on  standard  section  shapes  wars 
undertaken,  In  which  the  variations  of  outlet  flow 
snglo  and  total  pressure  loss  wore  measured  over  e 
range  of  Incidence.  Tests  of  this  kind  remained 
the  heels  for  blado  profile  selectloo  right  up  to 
the  I990's. 

The  advent  of  high  hypese  ratio  engine* 
necessitated  the  development  of  tho  single  stage 
transonic  fan,  which  was  very  difficult  both 
because  of  the  supersonic  relative  Inflow  on  tho 
outer  radii  and  because  of  tho  steaply  curved  flow 
path.  The  assumption  of  simple  radial  equilibrium 
and  the  use  of  traditional  subsonic  profiles  ware 
entirely  inadequate.’  Coplln*77'  presented  the 
history  of  transonic  fen  evolution  In  the  form  of 
Fig  5.  Before  1970,  thin  aerofoils  with  relat¬ 
ively  sharp  leading  edges  wars  Introduced,  but 
again  only  of  arbitrary  (double  and  later  multiple 
circular  arc)  shape.  A  major  advance  wss  mads 
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when  eh*  streamline  curvature  through- flow  me  hod 
was  Introduced.  This  enabled  blada  profile*  to  b« 
defined  along  stream  surface*.  As  can  be  seen  In 
Fig  3  this  i np roved  the  performance  of  M21 1-type 
fans  by  seme  IX  and  enabled  the  KB211-22B  to  be 
Introduced  at  a  satisfactorily  competitive  per¬ 
formance  level. 


FIGURE  5.  Transonic  Single  Stage  Fans 


It  was  not  until  the  last  few  years  that  the 
next  step  leproveeent  became  possible,  described 
In  Fig  5  as  "reduced  shock  loss",  as  a  result  of 
the  possibility  of  designing  high  speed  blade 
profiles  by  computational  methods  (rather  than 
using  circular  arcs).  The  principle  of  the 
"supercritical"  aerofoil  for  a  high  subsonic 
Inflow,  with  a  shape  tailored  to  provide  a  shock- 
free  slightly  supersonic  surf sc*  velocity  distri¬ 
bution  sad  a  delayed  diffusion,  had  been  evolved 
for  wings.  Sauer,  Carabsdlan,  and  Korn(78) 
published  a  computer  program  providing  the  corres¬ 
ponding  solution  for  a  cascade.  Coplln  reported 
the  result  of  applying  their  method  to  an  RB211 
fan  outlet  guide  vane.  The  number  of  vanes  was 
reduced  -  saving  weight  and  cost  -  and  the  effici¬ 
ency  Improved  by  \  to  IX  in  that  case. 

For  a  rotor  blada  section  with  supersonic 
relative  inflow,  no  purely  lavlscld  method  proved 
effective,  as  explained  earlier.  The  application 
of  Calvert's  method'32',  which  allowed  for  separa¬ 
ting  and  nearly  separating  boundary  layers  as  wall 
as  the  shock  pattern,  enabled  blading  with  lower 
shock  and  separation  losses  to  be  designed.  Two 
quite  different  fan  rotor  blades  (one  military  and 
one  elvll)  have  recently  been  designed  using 
Calvert's  aethod  and  tested  by  Kolls-ftoyce.  Both 
displayed  efficiency  Improvements  at  design  speed 
over  designs  undertaken  prior  to  the  availability 
of  the  new  computational  tool.  The  gains  ranged 
from  2  to  AX  and  were,  in  the  multistage  case, 
partly  associated  with  more  accurate  matching. 
Indeed,  this  is  a  good  emampls  of  how  a  well 
proven  computational  method  can  reduce  trlal-and- 
error;  only  one  build  was  needed.  This  first  time 
success  was  also  due  to  smother  computational 
Improvement ,  net  this  time  1m  the  fluid  dynamics 
area)  the  application  of  a  new  Rolls-toyce  method 
of  calculating  the  deformation  of  a  fan  blade 
under  centrifugal  and  aerodynamic  loading. 

The  successful  civil  rotor  desi|n  was  under¬ 
taken  by  a  proper  quasl-thrse-dlmensional  procsd ■ 
urs.  The  early  19t0's  designs  had  been  done  by 
first  calculating  the  through-flow  using  calculat¬ 


ing  planes  between  blade  rows  (not  within  them) 
and  then  blading  along  the  resulting  stream 
surfaces.  This  now  design  was  reached  by  Iterat¬ 
ing  between  the  blado-to-blade  calculation  (on 
stream  surfaces)  and  through-flow  calculations 
with  planes  within  the  blade  source,  implementing 
the  Ifu  approach  in  full.  Cinder^*3)  had  shown 
that  simply  Interpolating  conditions  within  the 
row  (knowing  those  between  the  rows)  could  load  to 
significant  errors,  as  the  radial  components  of 
blada  fores  are  than  unaccounted  for. 

CFO  has  played  a  less  significant  role,  so 
far,  in  ears  compressor  improvements.  The 
through-flow  is  first  calculated  by  a  streamline 
curvature  method,  using  pianos  between  the  blade 
rows  only.  With  meb  etralghter  annulus  wall 
lines  than  in  a  transonic  fan.  It  seams  lass 
necessary  to  consider  planes  within  the  rows:  and 
of  eourss  the  computational  grid  for  a  wit  is  tags 
compressor  might  become  too  large  with  the  extra 
planes. 

The  selection  of  blade  profiles,  for  many 
years  taken  from  "standard*  shapes,  has  recently 
followed  the  "supercritical"  route  previously 
described  for  fan  outlet  guide  vanes,  with  checks 
using  other  invlseid  methods*  Aerofoils  of  this 
type  have  consistently  shown  efficiency  Improve¬ 
ments  around  IX,  (equivalent  to  reducing  the 
actual  loss  by  the  order  of  10X)  and  reductions  In 
the  mother  of  blades  by  more  than  10X.  This 
reduction  in  blade  numbers  accounts  for  the 
efficiency  gain  end  economises  directly  In  coat 
and  weight.  The  new  shapes  could  not  have  been 
generated  by  any  simpler  method.  They  are  checked 
by  finite  difference  or  finite  element  methods. 

Centrlfuial  compressors 

The  early  OK  centrifugal  compressors  were 
designed  by  essentially  ona-dlmenslonal  methods. 
In  the  late  1960's  and  early  1970' s,  finite 
difference,  matrix,  and  streamline  curvature 
methods  bscsms  available  to  compute  impeller  vans 
surface  velocity  distributions  on  the  assumption 
of  unaeparatsd  flow.  RCTe(*3)  and  NELC**)  made 
their  libeller  design  schemes  evalleble  to  OK 
industry  In  the  lata  1970* s.  The  REL  methods  were 
applied  there  to  design  various  pmps  and  fans  for 
commercial  customers.  For  example,  a  quiet  cool¬ 
ing  fan  was  designed  in  1973  iriilch  bad  s  much 
batter  performance  than  its  competitors.  The  RCT8 
package  was  used  by  Comp  Air  to  design  a  most 
satisfactory  fan*7*',  and  was  adopted  and  tnpvod 
by  Noel  Penny  Turbinas  for  various  designs™0). 

The  application  of  scientific  design  princi¬ 
ples  to  centrifugal  pumps  can  have  a  startling 
effect  on  performance.  Fisher ^ 81 )  quotes  the  case 
of  an  automotive  water  pump  at  least  throe  times 
as  efficient  as  Its  prsdecsssor. 

In  the  aeronautical  field,  the  first  applica¬ 
tion  of  a  modern  through-flow  calculation  was 
worth  some  AX  In  efficiency.  A  mors  notable 
emampls  occurred  when  the  opportunity  can*  to 
tolls-Royce  to  redesign  tho  Dart  Impeller,  phlch 
had  originally  been  developed  from  pre-war  super¬ 
chargers  long  before  CFD  had  been  introduced.  The 
fuel  constBptlon  of  the  engine  was  reduced  by  a 
remarkable  g%X. 

It  Is  interesting  to  mote  that  all  the 
methods  actually  used,  so  far,  assumed  unaeparatsd 
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flow,  wheraes  latvt  —our— of  have  confirmed 
choc  most  Impeller  flows  appeared  co  have  local 
separations  M<(  of  the  time.  The  fully  thtee- 
dlatuloul  vlicow  KCbodo  now  hoc  owing  available 
should  eocor  for  oops roc loos.  Will  lc  become 
possible  then  to  design  even  bettor  impellers  end 
overcome  che  relatively  low  efficiency  of  eoec 
high  preeeure  redo  unite? 

Turbines 

Eerly  turbines  were  designed  assuming  uniform 
flow  from  hub  to  tip,  or  Inter  soon  form  of  rsdlsl 
equilibrium.  The  out  1st  gee  eagle  from  e  turbine 
blede  row  Is  felrly  well  approximated  by  cos*1 
( throat /pitch) ,  so  the  blade  profiles  were 
designed  on  a  drawing  board  by  narking  out  the 
throat  circle  (to  suit  che  required  outlet  angle) 
then  fitting  tewerel  circular  area  to  blend  into  e 
smooth  streamlined  shape.  Continuous  contraction 
of  the  passage  width  up  to  the  throat  wee  ensured. 

As  In  the  case  of  compressors,  the  eerly 
conformal  mapping  methods  were  not  ueod  to  design 
engine  blades,  chough  some  exact  solutions  gener¬ 
ated  by  these  methods  were  used  ee  test  cssss  for 
validating  approximate  maser leal  methods. 

C7D  methods  wars  first  applied  in  the  1960's, 
to  design  better  aerofoil  shapes.  Initially,  the 
incompressible  Mertensen  method^2*)  wee  used  at 
MCTE.  One  lesson  that  the  theory  llltsslnsted  wee 
that  surface  curvature  (not  elope)  le  the  geomet¬ 
rical  property  appearing  directly  in  the  squat lone 
determining  surface  velocity.  So  to  get  e  smooth 
surface  velocity  the  curvature  must  be  continuous. 
This  conflicts  with  the  Instinctive  feeling  that 
only  slope  needs  to  be  continuous.  So  profiles 
designed  using  contiguous  circular  arcs  did  not 
show  favourable  surf see  velocity  dlstrlbutloos. 
The  effect  of  using  the  Mar Cana an  approach  was  to 
Improve  efficiency.  In  one  case,  blading  designed 
by  NCTS  for  a  small  Industrial  turbine  manufact¬ 
urer  ebowed  e  7X  Improvement  over  e  traditional 
design. 

Rolls-Royce  adopted  the  Stenlts  PVD 
approach^**)  end  It  wee  widely  ueed  to  deelgn 
turbine  profiles.  The  method  designs  e  passage, 
using  e  linearised  compressible  flow  calculation; 
leading  end  trailing  edges  ere  added  afterwards. 
It  was  first  applied  to  the  Olympus  turbines,  and 
direct  comparison  between  old  and  new  designs 
(conducted  on  verlowe  configurations)  showed 
Immediate  galea  of  up  to  6X  le  efficiency*82) . 

In  the  1970' e,  tolls-Koyee  adopted  streamline 
curvature  methods  for  both  through- flow  sod  blede- 
to- blade  deelgn  end  analysis.  The  principle  of 
deelgn  wee  to  ecert  with  e  favourable  prosenre 
distribution  and  Iterate  manually  until  s  coolable 
•hepe  (that  Is,  an  sorofoll  with  sufficloat  thick¬ 
ness,  loading  and  trailing  edge  redll  end  trailing 
edge  wedge  angle)  was  reochod. 

In  the  1910*  s,  both  NOTE  end  Solle-torce 
turned  to  the  Denton  two-dimensional  program'**', 
which  enabled  locally  ewpersonlc  flows  (ee  no  nr 
ally  encountered  in  che  trailing  edge  region)  to 
be  calculated.  Unlike  e  compressor  passage,  e 
turbine  peeeege  le  not  greatly  effected  by  the 
relatively  chin  boundary  layers  on  the  blade 
surf sees  (oxeope  near  cho  annulus  walls).  So  an 
lnviscid  method  provides  a  good  basis  for  design. 
Thera  is,  unforconscely,  a  caveat  to  this.  In 


1983  Paige*7  ^designed  a  nossls  guide  vane  profile 
which  was  "bettor'*  chan  an  existing  profile  In 
that  the  pressure  distribution  looked  bettor  end 
calculation  of  the  surface  boundary  layers  led  to 
e  lower  loae  prediction.  When  tested  In  cascade, 
the  deelgn  velocity  distribution  wee  achieved  but 
the  overall  Lose  wee  nevertheless  higher  then  the 
old  design  because  the  bees  preeeure  wee  lower. 
The  beta  preeeure  cannot  be  predicted  by  an 
lnviscid  eat hod  end  the  change  wee  not  in  accord¬ 
ance  with  base  pressure  correlations.  This 
sxperlsncs  sounded  e  note  of  caution. 

In  the  lest  few  veers,  the  fully  three- 
dimensional  Denton  method*8*'  hee  been  employed  by 
1AI  end  by  Rolls-Eoye*  to  help  design  new 
eur>lnes.  The  proximity  of  an  annulus  wall  can 
significantly  change  the  aerofoil  surface  flow. 
The  availability  of  e  three-dimensional  method 
allows  the  deelgner  to  explore  changes  in  blede 
stacking  (the  way  aerofoils  st  each  radius  arc 
relatively  positioned)  and  changes  In  and  wall 
shape.  Morgen*8*)  has  described  the  effect  of 
changing  the  stacking  of  the  RB211  hp  nossle  guide 
vanes,  guided  by  C7D  methods,  which  increased 
efficiency  by  around  1Z. 

Another  method  used  In  recant  yeere  by  Rolle- 
Koyce  Is  the  Moore  three-dimensional  viscous 
program*8*).  This  is  formally  restricted  to  sub¬ 
sonic  compressible  flow,  but  It  can  be  used  to 
eesese  possible  secondary  flows.  A  particular 
example  In  which  the  lnviscid  and  viscous  three- 
dimensional  methods  were  used  to  guide  e  design 
change  was  the  RB211  IP  nossle  guide  vsne.  These 
venes  are  wounted  In  a  duct  of  rapidly  increasing 
radius,  which  cannot  adequately  be  catered  for  In 
e  quasi- 3D  method.  The  ngv  used  In  earlier 
engines  was  redesigned  to  remove  e  local  three- 
dimensional  flow  separation,  end  the  engine 
specific  fuel  consumption  Improved  by  soma  l\l. 

Turbines  ere  of  course  widely  used  out aide 
the  aircraft  Industry,  end  che  seme  CPD  methods 
are  avails bis,  for  example,  to  eteem  turbine 
designers.  It  has  recently  been  reported'83'1  that 
the  LP  turbine  In  one  of  the  three  Persona  500  MW 
seta  et  Dldcot  Power  Station  was  modified  by 
Persons  for  the  Central  Electricity  Generating 
Board  at  e  coat  of  £3M  to  Incorporate  e  new  lest 
stage  rotor  blede  designed  with  the  help  of  e 
Denton  30  program,  and  new  diaphragms*  The  unit 
efficiency  rose  some  3Z,  end  the  resulting  saving 
In  the  coal  bill  it  E1.7M  per  annum  on  that  one 
set  alone. 

At  the  other  end  of  the  else  range,  Connor 
end  Payne'88'  have  described  how  the  application 
of  FVD  methods  has  Increased  the  efficiency  of  e 
turbocharger  exlsl  carbine  by  over  3Z  et  design 
point,  increasing  to  over  10X  et  off-design. 

5  Comments  on  CPD  Klstorv 

Looking  beck  et  whet  has  happened,  e  number 
of  gonerel  commence  can  be  made: 

(1)  The  Invention  of  e  fundamentally  new  method 
Is  rsrt  and  la  done  not  by  e  teem  but  by  e  gifted 
Individual.  It  cannot  be  "scheduled"  by  e 
leeeereh  Manager;  ell  he  can  do  is  to  crests  the 
conditions  under  which  e  suitable  research  worker 
Is  attracted  to  eha  problem  end  equipped  with  the 
Clue  end  facilities  to  tackle  It. 


(2)  The  implementation  of  a  now  method  la  a  loot 
hard  Job  which  can  bo  greatly  helped  by  loaa 
gif tod  workers  than  cho  originator.  Vhon  It  comes 
to  converting  It  Into  a  "us or  friendly"  program 
capable  of  being  used  by  designers  who  do  not 
understand  the  mathematics,  there  la  a  great  deal 
of  work  to  be  done  by  a  com,  properly  planned  and 
professionally  managed. 

(3)  The  validation  of  a  new  prograw  agalnat 
experimental  results  la  vital  and  requires  a 
complicated  (and  probably  expensive)  experiment 
planned  and  executed  with  the  help  of  the  CFD 
analyst.  Especially  for  high  speed  turbo** 
aachlnery,  the  necessary  experimental  facilities 
can  only  be  found  at  national  Research  Establish- 
eencs  or  In  Industry.  Zt  must  be  the  function  of 
a  Research  Manager  to  plan  such  work. 

(4)  The  decision  to  commit  a  new  design  evolved 
by  CFO  methods,  first  to  an  experimental  demon¬ 
stration  and  later  to  a  production  anglne,  la  a 
Chief  Engineer’s  decision,  and  a  proposal  has  to 
be  Justified  to  him  by  unequivocal  validation 
achievements,  carefully  planned. 

(5)  Many  more  methods  have  been  evolved  than  have 
ever  been  put  to  good  use.  It  la  probably  not 
true  that  "only  the  fittest  survive".  Much  seems 
to  depend  on  the  accidents  of  history  s  which 
computer  wee  used,  how  eloquent  the  originator 
was,  which  organisation  he  happended  to  work  In. 
The  moat  marked  progress  occurred  when  two  or  more 
methods  were  actually  competing  scientifically  on 
the  same  problem.  Research  Managers  should 
encourage  such  competition  In  developing  new 
methods. 

(6)  The  visual  presentation  of  three-dimensional 
flow  patterns,  to  enable  the  research  worker  him¬ 
self  to  grasp  them,  and  later  to  enable  him  to 
explain  hls  results  to  others,  presancs  some 
difficulty.  Considerable  effort  has  needed  to  be 
devoted  to  computer  graphics. 

6  Future  Meed a  and  Prospects 

Fig  6,  showing  three  generations  of  combat 
aircraft  engine  scaled  to  the  same  thrust. 
Illustrates  the  overall  achievements  over  the  last 
twenty  years  In  engine  design  and  ehe  target  for 
the  immediate  future. 

Axial  compressors 

Despite  the  much  improved  ability  to  predict 
the  design  point  flow  field  and  performance  of  a 
transonic  compressor  stage,  described  earlier,  the 
reliable  prediction  of  off-design  performance  - 
and  especially  surge  prediction  -  remains  elusive. 
The  first  missing  element  la  the  accurate  predic¬ 
tion  of  the  total  preeaura  loaa  of  a  eharp- 
leadlng-edga  transonic  profile  at  off-deslgn 
conditions,  and  this  seems  achievable  by  the 
improvement  of  existing  methods.  A  ttaeh  more 
difficult  problem  -  requiring  e  fully  three- 
dimensional  viscous  method  -  la  the  prediction  of 
flow  separation  (possibly  leading  to  rotating 
stall  or  aurga)  in  the  end  wall  regions,  including 
tip  clearance  effects.  Another  reason  for  needing 
a  3D  viscous  method  is  to  predict  the  radial 
migration  of  aerofoil  boundary  layer  fluid  which 
effectively  "transfers"  loaa  from  one  radius  to 
another,  and  may  on  occasions  accumulate  low 
•nargy  fluid  around  part-span  shrouds. 


SPEY  1969  21  Stapes 


RB199  1981  16  Stages 


FUTURE  Late  1990  a  9  Stages 


FIGURE  6.  Combat  Engines  -  Scaled  to  the 
Same  Thrust 


The  trend  of  compressor  detign  has  been  not 
so  much  towards  Increasing  efficiency  as  maintain¬ 
ing  efficiency  at  progressively  incraaslng 
praasure  ratio  par  stage.  Pig  6  shows  how  this 
trand  has  reduced  the  elze  and  weight  of  military 
engines  and  la  continuing  to  do  so;  the  three 
engine  drawings  illustrated  have  been  scaled  to 
the  same  thrust.  The  principal  element  In 
Increasing  atage  pressure  ratio  la  of  course 
Increasing  blade  apeed  as  Improved  materials  ara 
developed;  but  the  Inevitable  consequence  le  more 
and  more  supersonic  flows.  If  a  military  engine 
la  ever  to  have  a  single  atage  fen  to  achieve  Its 
pressure  ratio  of  3  or  4,  it  would  be  a  fully 
■uperaonlc  compressor.  A  great  deal  of  effort  waa 
devoted  to  abortive  attempts  to  produce  efficient 
supersonic  compressors  in  the  1950's  and  1960's; 
could  the  CFD  tools  now  being  developed  enable 
much  more  successful  attempts  to  be  mede  in  the 
1990'a? 

Turbomachine  flows  are  assumed  to  be  steady 
In  moat  analyses  and  designs.  Unsteady  calcula¬ 
tions  have  been  concerned  primarily  with  flutter. 
It  la  not  yet  clear  whather  the  development  of 
unsteady  flow  computations  (as  computers  become 
fast  enough)  will  reveal  e  need  to  alter  design 
philosophies  to  Improve  performance.  Reliable  CFD 
predictions  of  flutter  and  noise  appear  very 
difficult  end  distant  targets  still. 

Turning  to  multistage  core  compressors,  there 
la  again  a  difference  between  the  priority  targets 
for  civil  engines  (efficiency)  and  military 
anginas  (compactness  and  low  weight),  but  the 
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technology  required  to  achieve  then  ie  common. 
The  flow  in  the  end  wall  regions,  including  the 
tip  clearance  flow,  has  a  doainane  Influence  on 
the  annulus  vail  boundary  layer  blockage  and  hence 
stage  Matching.  It  also  plays  a  large  part  in 
generating  total  pressure  losses,  and  perhaps  in 
surge  initiation.  Three-dimensional  viscous 
progress  are  needed  to  give  the  designer  a  better 
insight  into  the  phenomena,  and  especially  to 
guide  the  evolution  of  "end  bends". 

The  iaportance  of  nixing  within  a  Multistage 
compressor,  by  both  secondary  flow  and  turbulent 
eddies,  has  been  explained,  and  further  research 
in  this  area  should  contribute  substantially  to 
design  techniques. 

In  addlt.  on  to  improved  performance  in  the 
final  version,  the  application  of  CFD  will  signi¬ 
ficantly  reduce  the  number  of  trial  builds 
required  in  development,  so  saving  time  and  money. 

Finally,  the  response  of  compressors  to  a 
distorted  inlet  flow  (typically  due  to  combat 
•anoeuvres)  has  been  extensively  studied  experi¬ 
mentally  but  to  date  only  two-dimensional 
theories'®7'  (allowing  circumferential  and  axial 
flow  variation  but  no  radial  variation)  have  been 
employed  successfully  for  analysis.  A  three- 
dimensional  theory  (allowing  radial  variation)  Is 
clearly  essential  for  a  low  hub/tip  ratio  trans¬ 
onic  fan,  and  possibly  also  for  a  core  compressor. 

Attempts  to  date'52,80)  h4Ve  produced  methods  too 
restricted  in  scope  or  too  difficult  to  apply.  A 
solution  of  this  problem  appears  possible  using 
modern  msBcrlcal  methods  and  should  be  attempted. 

Centrifugal  compressors 

In  the~  1960's,  extensive  research  was 
devoted,  especially  in  Canada,  to  high  pressure 
ratio  units,  which  were  successfully  developed 
only  at  the  expense  of  lower  efficiency.  The 
losses  appeared  Insvl cable  because  of  the  high 
supersonic  Inflow  to  the  narrow  diffuser  ring. 
Most  small  aeronautical  gas  turbines  have  there¬ 
fore  chosen  to  use  several  axial  stages  followed 
by  a  lower  pressure  ratio,  more  efficient,  centri¬ 
fugal  stage.  The  other  problem  of  high  pressure 
ratio  centrifugal  stages  was  of  course  stressing, 
but  improvements  In  materials  will  presumably 
continue.  So  there  appear  to  be  two  llnaa  of 
advance  for  centrifugal  compreseore,  both  heavily 
dependent  on  CFD  lmproveaente. 

The  flret  le  the  low  pressure  ratio  unit  - 
around  3  -  when  the  flow  reechea  the  diffuser 
subeonlcally.  As  already  aentlonad,  the  Impeller 
flow  le  probably  asperated,  end  it  seems  e  reason¬ 
able  target  to  evolve  unseparated  designs  end 
hence  increase  overall  efficiency  to  axial  com- 
praaaor  lavela.  A  centrifugal  stage  could  then 
become  attractive  even  for  e  large  civil  engine. 
The  diffuser  Introduces  eonaldarabla  loss  because 
the  wetted  surface  area  la  larga,  and  thara 
appears  to  be  scops  for  tha  application  of  3D 
viscous  codas  to  improva  dlffuaars. 

Tha  ascend  lino  of  sdvsncs  could  b#  s  return 
to  higher  pressure  ratio  unite  (10?)  to  reduce 
engine  weight  and  coat,  but  reducing  the  effic¬ 
iency  lose  by  tailoring  the  shock  patterns  with 
the  help  of  3D  codes. 


Axial  turbines 

Al though  the  equations  ere  the  sene,  there 
are  significant  dlfferencea  between  compressors 
and  turblnas  for  tha  CFD  analyst.  The  first  is 
that  the  turbine  blade  surface  boundary  layers  are 
small  fractions  of  the  passage  width,  so  that 
lnviscld  methods  give  e  good  guide  to  optimum 
shape  end  predict  the  flow  well,  until  near  the 
trailing  adge.  The  second  difference  la  that  the 
blade  camber  and  hanca  the  secondary  flows  are 
very  much  grsstar  than  In  compressors,  and  era 
accentuated  by  the  tip  clearance  rather  then 
reduced.  A  third  difference  Is  that  hast  transfer 
le  a  key  element  in  turbine  design. 

Most  aeronautical  turbine  nossle  guide  vans 
rows  and  many  rotor  rows  have  sonic  or  supersonic 
relative  outflow,  with  s  shock  structure  Impinging 
oo  the  boundary  layer  Just  ahead  of  tbs  trailing 
edge.  The  resulting  lmmbds-shock  controls  the 
best  pressure.  Currently,  the  bass  pressure  le 
predicted  empirically  and  It  seams  possible  that, 
oo  a  two-dimensional  basis,  a  viscous- lnviscld 
interact ivs  method  or  a  fully  viscous  method  might 
be  able  to  predict  It  theoretically.  lecause 
turbulence  modelling  of  separated  regions  le  still 
difficult,  the  interactive  method  seems  more 
likely  to  succeed  in  the  short  term.  However, 
there  are  two  serious  comp  lice  Clone;  one  is  the 
n«md  to  model  tome  form  of  radial  equilibrium  In 
the  nearly  stagnant  base  region,  and  the  other  is 
the  effect  of  cooling  sir  discharge  into  the  bees 
region.  Both  phenomena  appear  to  offer  scope  for 
Improved  overall  efficiency  if  they  could  be  well 
enough  understood. 

Considerable  quantities  of  experimental  date 
have  been  amassed  on  end  well  end  secondary  flow. 
The  secondary  flow  tends  to  etrip  the  incoming 
well  boundary  layer  fluid  off  the  wall  and 
discharge  it  into  the  mainstream  via  the  blade 
suction  surface  trailing  edge.  A  new  wall 
boundary  layer  etarta.  The  ability  to  predict 
this  type  of  flow  depends  upon  the  development  of 
e  3D  viscous  code  able  to  ester  for  the  corners 
between  well  and  binds,  and  In  dus  course  also  for 
tip  clssrsncs.  The  xero  clearance  ceee  seems 
likely  to  be  solved  quite  soon.  It  Is  gsnarally 
accepted  that  end  wall  effects  of  this  kind 
account  for  around  half  the  total  prsssurs  loss  in 
s  turblns;  sursly  that  loss  can  be  reduced  by 
sclentifleally-choeen  end  wall  profiling  or  "and 
bending"  of  the  blades,  when  those  CFD  tools 
become  available.  For  a  multistage  turbine,  the 
type  of  mixing  analysis  described  for  axial 
compressors  le  surely  also  needed. 

The  prediction  of  external  heat  transfer  Is 
largely  e  question  of  boundary  layer  prediction, 
which  becomes  particularly  difficult  In  the 
preeence  of  film  cooling  or  GQrtler  vortices. 
However,  It  haa  been  demonstrated  that  the  passage 
of  upstream  wakes  (or  even  of  downstream  blades) 
has  a  aalor  effact  on  the  boundary  layers (M,M). 
Hobson^*!)  showed  the  boundary  lsysr  switching 
from  laminar  to  turbulant  and  back  again  as  wakes 
passed.  It  is  clear  that  an  analysis  of  this 
situation  raqulrss  an  unsteady  CFD  model.  At 
least  on  a  two-dimensional  viscous  basis,  this 
should  ba  already  possible. 

As  for  compressors,  CFD  will  not  only  Improve 
turbine  performance  but  also  reduce  development 
time  end  cost. 
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Radial  tgrblMi 

tad  la  X  ieflow  turblnea  have  not  bean  used  In 
OK  eero  engines,  not  on  account  of  inefficiency 
but  because  of  tba  mechanical  design  problem  of 
asking  a  cooled  high  speed  rotor  of  adequate 
strength  and  life,  sod  not  too  high  a  rotational 
Inertia.  Aa  a  result,  they  have  not  received  the 
attention  of  many  CTO  specialists  In  the  UK. 
There  must  be  scope  for  applying  the  methods 
developed  for  centrifugal  compressors,  without  the 
worries  about  large  local  flow  separations*  Sons 
work  of  this  kind  has  been  done  for  application  to 
large  turbochargers^6). 

The  possibility  of  advanced  eeraalc 
materials,  not  needing  cooling,  could  promote  the 
radial  turbine  as  an  option  for  small  aeronautical 
applications. 


7  Conclusions 

Computational  fluid  dynamics  methods  specific 
to  turbomachinery  have  been  developed  since  1940 
within  the  UK.  They  have  usually  employed  the 
same  mathematical  methods  used  In  other  branches 
of  CTO,  aa  and  when  computers  became  large  and 
feat  enough.  Initially,  the  models  used  were  too 
simplified  (incompressible  Invlscld  planar  or 
cylindrical  flow)  to  be  realistic  end  too  labor¬ 
ious  to  apply  end  were  not  used  in  practice.  The 
practicability  of  achieving  realistic  design  and 
analysis  results  using  the  quasi- three-dimensional 
(Wu)  approach  meant  that  CTO  could  however  be 
usefully  employed  aa  early  aa  the  1960's  although 
computers  than  could  only  calculate  two- 
dimensional  flow  fields.  By  the  1980's,  three- 
dimensional  Invlscld  flow  fields  could  b*  computed 
and  the  present  decade  Is  seeing  great  advances 
Including  viscous  calculations. 

The  UK  aeroengine  firm  (Kolle-toyce)  has  been 
quick  to  adopt  CTO  methods  sod  to  employ  them  to 
design  and  develop  better  compressors  end 
turbines.  Just  as  other  aeroengine  manufacturer* 
have  don*.  The  first  bay  requirement  for  this  to 
happen  la  that  the  firm  should  have  a  "cor*  CTO 
teem"  large  enough  not  only  to  develop  its  own  new 
methods  sometimes,  but  essentially  to  take  s 
chosen  method  from  any  source  and  develop  It  Into 
a  proven  working  tool  capable  of  use  by  s  designer 
who  does  not  understand  the  mathematical  or 
programming  details.  The  second  requirement  is 
that  rasesrch  managers  la  the  firm  or  associated 
Research  Establishment*  must  organise  a  systematic 
methods  validation  programme,  capable  of  convinc¬ 
ing  the  moat  hard-headed  Chief  Engineer  to  coalt 
a  CTO  design  to  his  angina. 

Examination  of  ehe  history  of  UK  turbo- 
machinery  technology  since  I960  shows  that  some  of 
the  major  advances  In  product  quality  were  made  aa 
a  direct  result  of  the  application  of  CTO,  which 
arguably  could  not  have  been  achieved  without  It. 
Typically,  I960' a  methods  Improved  on  1950 'a 
nethods  by  soma  5  to  10Z  in  efficiency,  that  Is, 
reducing  losees  by  JO  to  50Z.  The  subsequent 
improvements  have  Increased  the  stage  loading 
levels  at  which  high  efficiency  can  be  maintained. 
Thre^-dlmanalonal  flow  field  tailoring  -  only 
imperfectly  understood  and  not  yet  predictable  in 
CTO  tarns  -  has  proved  generally  worth  another  1 
to  2Z. 


The  future  trend  In  CTO  le  Inevitably  towards 
3D  viscous  flow,  sod  unsteady  effects.  As 
computers  improve,  these  new  methods  mill  make 
available  to  the  turbomacMnery  designer  on  a  more 
rational  basis  *  vide  range  of  options  primarily 
In  the  and  wall  regions:  wall  profiling,  and 
bands,  varying  stacks,  wlnglete,  tip  treatment, 
easing  treatment.  Thar*  is  obvious  scope  for 
reducing  end  wall  region  losses  which  amount  to 
around  half  th*  total  losses.  The  possibility  of 
significant  savings  in  the  number  of  blades  and 
stages  -  and  banc*  In  weight  and  coat  -  will  art  e 
from  a  batter  understanding  of  how  to  control 
supersonic  flows*  Finally,  the  escalating  cost  of 
engine  development  may  be  reduced  considerably  by 
getting  ch*  aerodynamics  ''right  first  time"; 
visible  progress  has  already  been  made  in  this 
direction. 

Examples  of  successful  application  of  CTO 
methods  to  oon-aeronautical  turbomachines  have 
also  been  given.  There  are  two  industries  -  aero¬ 
space  and  power  generation  -  with  a  major  economic 
justification  for  performance  Improvement  through 
undertaking  CFO  research*  Th*  methods  they  gener¬ 
ate  will  continue  to  be  adopted  by  other  turbo- 
machinery  Industries. 

It  la  concluded  that  the  further  advance  of 
computational  fluid  dynamics  for  turbomachinery 
should  be  vigorously  encouraged  because  it  has  a 
najor  role  to  play  In  advancing  technological 
standards  and  reducing  th#  time  and  coat  of  devel¬ 
opment. 
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